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Abstract 

Azimuthal distributions of charged particles and intermediate mass frag- 
ments emitted in Au+Au collisions at 600^ MeV have been measured using 
the FOPI facility at GSI-Darmstadt. Data show a strong increase of the in- 
plane azimuthal anisotropy ratio with the charge of the detected fragment. 
Intermediate mass fragments are found to exhibit a strong momentum-space 
alignment with respect of the reaction plane. The experimental results are 
presented as a function of the polar center-of-mass angle and over a broad 
range of impact parameters. They are compared to the predictions of the 
Isospin Quantum Molecular Dynamics model using three different parametri- 
sations of the equation of state. We show that such highly accurate data 
provide stringent test for microscopic transport models and can potentially 
constrain separately the stiffness of the nuclear equation of state and the 
momentum dependence of the nuclear interaction. 

Keywords : Azimuthal distributions, nuclear equation of state, transport models, incom- 
pressibility, momentum dependent interaction. 

PACS numbers : 25.75.Ld 
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I. INTRODUCTION 



One of the major objectives in the investigation of high energy heavy ion colhsions 
(-Eiab > lOOAMeV) is to provide information on the equation of state (EoS) of nuclear 
matter under extreme conditions of density and temperature. In this respect, collective 
flow effects of various ejectiles emerging from the reaction are of great interest since they are 
known as a sensitive probe of compression of nuclear matter [T]] . These collective phenomena 
were predicted by nuclear hydrodynamical models and experimentally identified in the 
eighties with the advent of the first generation of Air detectors (Plastic Ball |^ and Streamer 
Chamber § at LBL-BEVALAC and the DIOGENE detector § at SATURNE) capable of 
full event by event characterisation. Different forms of nuclear matter flow were observed : 
i) the directed sideward flow which appears as a sideward deflection of the emitted particles 
in the reaction plane ||^,|^ ii) the squeeze-out effect which manifests itself as an enhanced 
emission out of the reaction plane and iii) the so-called radial flow which has been 
more recently evidenced by the FOPI collaboration as an azimuthally symmetric expansion 
in highly central collisions . 



Despite the extensive experimental studies |[TT| and the substantial theoretical 
progress p!2|-pO[] achieved so far, an unambiguous determination of the nuclear EoS from 
heavy ion collision studies is not yet at hand to date. Indeed there remain still several 
difficulties in addressing this intriguing question : 

• The extraction of information from experimental observables demands a thorough un- 
derstanding of the reaction dynamics which can be possible only through prerequisite 
tests of nuclear transport models. 

• The experimentally reconstructable observables are not uniquely sensitive to the in- 
compressibility coefficient (K) characterising the nuclear EoS. Other effects such as the 
momentum dependence of the nuclear interaction (MDI) and the in-medium modifica- 



tion of the nucleon-nucleon cross section ((Tnn) play also an important role |T2|,|T3|,|T^,pT . 
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It is therefore crucial to disentangle the influence of the EoS from the one related to 
the other physical effects above mentioned. 

There are also several input parameters needed for the initialisation of transport 
codes and somewhat arbitrarly fixed in the calculations, which might also in- 
fluence significantly the theoretical predictions. 

The sensitivities of the experimental flow observables to the EoS (as well as MDI and 



in-medium cTnn) are generally quite low [p!3ip3| , p4| , which requires highly accurate data. 



This implies in particular a good accuracy on the reaction plane reconstruction and 
low detector biases and inefficiences. 



Available transport models underpredict composite particles yields [p5H 28|. Thus, a 
meaningful comparison between experimental flow observables and those predicted by 
the theory should be made on the basis of coalescence invariant quantities including the 
contributions of all emitted particles. This requires in particular the inclusion of inter- 
mediate mass fragments (IMF, Z > 3) which are known to carry larger amount of flow 
than do lighter particles p5|-|52| and which constitutes a non negligeable fraction of the 
total mass of the system even at incident energies of a few hundred A MeV [P^,P^-pi 



Several experimental observables, corresponding to different projections of the triple 
differential cross section, have been introduced in order to characterise quantitatively the 
strength of the collective flow in heavy ion collisions ||8| j35|j36[] . One of these observables is 
the so-called "in-plane azimuthal anisotropy ratio" (defined in section IV) which can be 
obtained from measurements of azimuthal distributions relative to the reaction plane. The- 
oretical calculations in the framework of the Boltzmann-Uehling-Uhlenbeck (BUU) 



model have shown the relevance of this observable as a sensitive quantity to pin down 
the nuclear EoS. Comparisons between data and BUU calculations were carried out for neu- 
trons pU|p][ , but in this case the sensitivity to variations of the equation of state was found 



to be rather low. Recently, it was reported |^ that IMF, because of their strong alignment 



along the flow direction p9|, correlate much better with the value of K used in the theory. 

We report in this paper the first detailed experimental results on the in-plane azimuthal 
anisotropy ratio of light charged particles and IMF. The data were obtained from azimuthal 
distribution measurements carried out for the Au+Au reaction at several bombarding ener- 
gies from 100 to 800A MeV, using the FOPI facility at the SIS (GSI-Darmstadt) accelerator. 
Here we focus on the high incident energy of 600A MeV where : i) the participant and the 
spectator products are better separated, ii) compressional and MDI effects are expected to 
be large |[T2|Jl7| and iii) high statistics dynamical model calculations are available We 
find that the in-plane anisotropy ratio increases with the charge of the detected fragment. 
We show also that IMF exhibit a very pronounced alignment along the flow direction. The 
dependencies of the anisotropy ratio on the polar angle in the center-of-mass (cm.) frame 
and the collision centrality are presented. They are compared to the predictions of the 
Isospin Quantum Molecular Dynamics (IQMD) |^ model for different parametrisations 
of the equation of state. Simulated events are filtered through the FOPI acceptance and 
sorted in accordance to the adopted experimental procedure, which allows direct and re- 
alistic comparisons between data and model predictions. We show that the availability of 
such experimental data provide a stringent test for dynamical transport models. Based on 
FOPI data it has been shown that the shape of the fire-ball in central collisions depends 



sensitively on a^n [111- Here, we emphasize the possibility to disentangle compressional and 
MDI effects on the basis of detailed azimuthal distributions in non central collisions for the 
participant as well as for the spectator source. 



II. EXPERIMENTAL SETUP 

The data were obtained in a complete experimental study of the Au+Au system at 6 
beam energies going from 100 to 800A MeV at the SIS (GSI-Darmstadt) accelerator using 
the FOPI/phase-I facility. As already noted, the present work is restricted to the incident 
energy of 600^4 MeV. At this energy, about 10^ events were collected under the central trigger 
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condition [^. The latter was defined by adjusting the charged particle multiplicity to a 



value which corresponds to impact parameters less than 9/m in a clean-cut geometrical 
model. Most of the data which will be presented in the following were obtained from the 
analysis of these central events. Only those related to the centrality dependence will include 
events taken under the minimum bias condition {i.e. events where at least 4 charged 
particles were detected in the angular range 7° < 0iab < 30°). This large amount of available 
events allowed us to extract high statistics data. 



The FOPI detector has been described in details in reference Here, we recall briefly 
its different components and main features. At the incident energy considered in this paper, 
the setup consisted of a highly segmented Forward Wall of plastic scintillators, covering in 
full azimuth the laboratory polar angles (Oiab) from 1.2° to 30°, complemented by a shell 
of thin plastic scintillator paddles mounted in front of it and subtending the Giab = 1-2° - 
7° angular range. This device allowed us, event by event, to identify the nuclear charge and 
measure the vector velocities of most of the light charged particles and IMF (up to Z = 12) 
emitted in the forward cm. hemisphere. Its high granularity (about a factor of three higher 
than the one of the earlier Plastic Ball detector 0) allowed high multiplicity events to be 
measured with a low multi-hit rate. The apparatus provides detailed and highly accurate 
measurements of triple differential cross sections for all charged particles with a very good 
accuracy on the reaction plane determination which is of crucial importance in the present 
work. Figure |l] illustrates the acceptance of the Phase-I configuration of the FOPI detector 
in the {pi^\y^^^) plane, where p[^^ is the transverse momentum per nucleon normalized to 
the cm. projectile momentum per nucleon (= 528.6 MeV/c) and y^^^ denotes the rapidity 
normalized to the projectile rapidity in the cm. frame. Solid curves represent the polar 
geometrical limits of the FOPI apparatus in the laboratory frame (Oiab = 1-2° and 30°). 
The detection thresholds are indicated by the dotted curves in the case of Z = 3 fragments. 
The dashed curves mark different cm. polar angles (Oc.m.) separated by 10 degrees. As 
it can be seen in Fig. |I|, the apparatus allows undistorded measurements (free of detector 
biases) to be performed over a broad region of the forward cm. hemisphere. 



III. EVENT CHARACTERISATION 



The study of azimuthal distributions requires a good event characterisation both in 
impact parameter and azimuthal orientation of the reaction plane. This is also of crucial 
importance when comparing the experimental data to theoretical predictions. In this respect, 
it is interesting to point out that a heavy system such as Au+Au has the advantage to 
offer a better determination of the collision geometry, both in centrality and reaction plane 
orientation, than lighter systems. In this section we will describe the procedure used in 
in order to characterize the measured events and illustrate the performances of the FOPI 
apparatus. 

To sort out the events according to their degree of centrality, we have employed the stan- 
dard method based on the correlation between the multiplicity of the emitted particles and 
the impact parameter. The multiplicity distribution of charged particles measured in the 
angular range 7° < 6iab < 30° exhibits a typical plateau for intermediate values and falls off 



for the highest multiplicities p6| , p3| . This distribution was divided into five intervals, in ac- 
cordance to the procedure introduced in previous works p3| , ^ . The highest multiplicity bin 
(named PM5) starts at half the plateau value corresponding to ^ 2.5% of the total reaction 
cross section and a maximum impact parameter of ~ 2.5/m in a sharp cut-off approxima- 
tion. The remaining multiplicity range was divided into four equally spaced intervals named 
PMl to PM4. The limits of these event classes are given in Tab. |. In the following only 
events belonging to the PM2-PM5 bins are considered. An additional condition requiring 
the total detected charge to be larger than 30 was imposed in the analysis, in order to reduce 
the contribution of the background contamination from non target collisions. Within this 
condition the background contamination was estimated, by comparing measurements with 
and without target, to be less than 5% in the PM2 multiplicity range and negligible for 
higher multiplicity events. The mean impact parameters and the associated r.m.s. devia- 
tions listed in Tab. |T| were extracted from IQMD simulations, where the theoretical events 



7 



were passed through the detector acceptance and sorted in accordance to the above outhned 
procedure used for the data. By doing so, one can compare directly the experimental ob- 
servables to the the theoretical ones, even if the calculations do not reproduce the measured 
multiplicity distribution, which is the case of the IQMD model. The latter, as we will see 
in section V, underpredicts IMF multiplicities and overpredicts the overall charged particle 
multiplicities. This can be seen in Tab. |, where the lower limit of the PM5 bin for IQMD 
events is much larger (by 14 units) than the experimental one. However, when scaled to 
the lower limit of the PM5 bin, the theoretical multiplicity distribution is found to be in 
a fairly good agreement with the data in the range PM2-PM5. As shown in Tab. |, the 
charged particle multiplicity criterium covers a rather wide impact parameter range from 
< b >~ 3fm to < b >~ 9/m. 

The reconstruction of the azimuth of the reaction plane have been performed according 
to the transverse momentum analysis |3^. In the framework of this method, the reaction 



plane is determined, for each particle /i in a given event, by the vector Q calculated from 
the transverse momenta of all detected particles except the particle in order to remove 
autocorrelation effects : 

M 

Q = J2 ^uPt 



I/=l 



M is the multiplicity of the event and ujy = \ \i y^^^ > S, —1 if y^^^ < 5 and otherwise. 
6, choosen equal to 0.5 in the present work, is a parameter introduced to remove mid- 
rapidity particles which are less correlated with the reaction plane. In order to estimate 
the accuracy on the reaction plane determination, due to finite number of particle effects 
and detector biases, we have used the method described in reference which consists in 
randomly subdividing each event into two and calculating on average the half difference in 
azimuth (A$r) between the reaction planes extracted from the two sub-events. A$r gives 
an estimate of the dispersion of the reconstructed reaction plane with respect to the true 
one [^. The results are displayed in Tab. | in terms of the standard deviation width o"(A$r) 
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extracted from a gaussian fit to the A$r distributions. As can be seen, tlie reaction plane 
is, in all cases, rather well estimated, with a precision which varies typically from ~ 24° to 
~ 34° depending on the event multiplicity. The data presented in this paper are corrected 
for these uncertainties on the reaction plane determination (see below). 

IV. EXPERIMENTAL RESULTS 

Azimuthal distributions are generally used to study the emission pattern of particles 
in the participant region {i.e. at mid-rapidity) of the collision. Instead of restricting the 
analysis to a narrow rapidity window centered around the cm. rapidity, we have explored 
here the whole experimentally covered phase space (see Fig. |l]). To select particles emitted 
in different regions of the phase space the analysis was carried out by imposing cm. polar 
angle (6c.m.) gates to the data, as suggested in previous works [^^. The use of this 
Oc.m. binning was found to be more suited than using rapidity cuts because of the typical 
FOPI/phase-I acceptance (Fig. |1|) which is limited to laboratory polar angles lower than 30°. 
Indeed by doing so, one can extract an accurate quantitative information from the measured 
triple differential cross sections over a broad region of the forward cm. hemisphere where, 
as we will see in the following, the data are only very weakly affected by the detector cuts. 

Typical examples of azimuthal distributions around the beam axis Y{^) = dN/d^ (where 
$ is the azimuth of the emitted particle with respect to the azimuthal orientation of the 
reaction plane) of charged particles are displayed in Fig. ^ for the PM4 event class. These 
distributions were obtained for proton-like particles, i.e. by including all detected particles 
each being weighted with its charge. Different panels in Fig. ^ show results for different 
cm. polar angle gates each 10 degrees wide except for the one around 90° which was 
taken from 80° to 100°. With increasing Oc.m., a drastic change in the azimuthal emission 
pattern is observed. At forward polar angles, the distribution exhibits a strongly enhanced 
in-plane emission along the sideward flow direction ($ = 0°), while at large polar Oc.m. 
angles (approaching Oc.m. = 90°), one observes a clear preferential emission in the out-of- 
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plane ($ = 90° and $ = 270°) direction corresponding to the squeeze-out effect. Thus, the 
shape of these azimuthal distributions reflects the strength of the collective motion for both 



the directed in-plane flow and the out-of-plane emission. In a previous paper we have 
reported the results concerning the preferential out-of-plane emission. Here we focus on the 
in-plane anisotropy component. 

In order to examine the data in a quantitative manner, we have extracted the magnitude 
of the in-plane anisotropy R for each Oc.m. bin as : 



^ _ \^ ) l0°<'I'<45° \^ ) l315°<'I'<360° ^-|^'^ 





0°<<I><45° 




315°<<I><360° 




135°<'I'<225° 



This deflnition of the in-plane anisotropy signal differs from the one used in earlier 
works pO| , ^ , namely the ratio |$^qo / Y{^) |$=i8o°- '^^^^ dictated by statistics 



limitations for the theoretical IQMD calculations. 

Fluctuations of the reconstructed reaction planes with respect to the true ones cause an 
attenuation of the experimentally observed anisotropics. To take into account this effect, 
we have proceeded in the following way 



• We have flrst fltted a function of the form oq + ai cos($) + 02 cos(2$) to the measured 

distributions. The resulting flts are shown by the curves in Fig. ^ 

• Then following we have determined the corrected parameters : 

a[j = ao , a; = ai / (< cos(A$r) >) , ^ = aa / (2 < cos2(A$r) > -1). 

• Finally, we have extracted the corrected anisotropy ratios R by integrating in equation 
(1) the corrected fitting function. 

The values of cos(A$r) and cos^(A$r) were determined event wise and then averaged 
for each event class. The results are shown in Tab. |T| in the case of the 40° < Oc.m. < 50° 
gate for different centrality bins. As can be seen, the corrections become more important 
for peripheral events (PM2 bin) because of the relatively large dispersion in A<I>r. 
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Figure ^ shows the in-plane anisotropy ratio R extracted from the data according to 
the definition in equation (1) and corrected for fiuctuations of the reconstructed reaction 
plane, as a function of the polar angle 6c.m.- The results were preselected by the charged 
particle multiplicity cut PM4. They are presented for different types of particles {Z = 1 
to 4). The observed trend is characterized by a well defined maximum located at forward 
cm. polar angles (Oc.m. — 15°). The presence of this maximum is related to the bounce- 
off of the projectile remnants. For polar angles approaching Oc.m. = 90°, the anisotropy 
ratio tends to R =1 as expected from symmetry reasons. Here, the azimuthal distribution 
exhibits a completely different pattern with an enhanced out-of-plane emission (Fig. whose 
magnitude is commonly characterized by the squeeze-out ratio P,^. It should be pointed 
out that the overall trend, observed in Fig. |, is qualitatively consistent with earlier data ||5^ 
obtained at much lower beam energy. 

As can be seen in Fig. ^ with increasing fragment charge one observes larger azimuthal 
asymmetries with a more pronounced maximum whose location is shifted towards smaller 
polar angles. These features are more quantitatively illustrated in Fig. ^, where the maxi- 
mum (top-panel) of the correlation R versus Oc.m. and its location O™^^ (mid-panel) 
are plotted as a function of the charge of the detected fragment, i?™*^^ and O™^^ were ex- 
tracted by fitting a gaussian to the measured correlation between R and logiQ{Qc.m.)- The 
resulting fits are shown, in the angular range where the fitting procedure was applied, by 
the dotted curves in Fig. |^. We have also displayed in Fig. ^ (lower-panel) the mean frag- 
ment multiplicities in the forward cm. hemisphere (y^^^ > 0). It is worth to note that 
IMF are still observed with quite large multiplicities at the relatively high incident energy 
of 600AMeV explored in the present work. We found |51|] that IMF carry about 12% of 
the total detected charge in the case of the PM4 event class. It should be noted that for 
this class of events, a substantial fraction of IMF originates from the decay of the spectator 
source. When the latter source is selected the fraction of IMF remains nearly the same in 
the whole beam energy range from 250 to SOOAMeV what is in agreement with the 
universal behaviour of the spectator decay reported by the ALADIN collaboration The 
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observed dependence between and Z (top-panel in Fig. ^) shows a roughly linear rise 
of the maximum anisotropy ratio with increasing fragment charge at least up to Z = 4. 
Above Z = 4, statistics in the data was not sufficient to extract accurate R values, because 
of the very large anisotropies and the low multiplicities. On the other hand, the mid-panel 
of Fig. ^ shows that with increasing Z the location O™^^ of the maximum anisotropy moves 
towards lower values and seems to saturate. 

The observed dependencies of B^^^ and Q^^^ with the fragment charge, reflect the 
strong alignment of IMF along the flow direction, an effect which was qualitatively observed 



in the earlier Plastic Ball experiments ||29| for the same Au-|-Au system at a beam energy 
of 200A MeV. They are consistent with the increase of the sideward flow with the fragment 
size which was also reported in reference and recently confirmed in more quantitative 



studies by several groups ||30|- p2| , pT|^ . The strong alignment of IMF could be qualitatively 



understood [^,0 from a simple picture where the pure collective motion is superimposed 
upon the random thermal one : light particles being more sensitive to thermal fluctuations 
which tend to wash out their flow, while IMF being less affected by the undirected thermal 
motion are expected to be much more aligned along the flow direction. 

V. COMPARISON TO IQMD MODEL PREDICTIONS 

In this section we compare the experimental results to the predictions of the QMD 
model [p!3tp!7| , ^p5| using the so-called IQMD version Besides the test of the IQMD 



model, our main aim here is to examine the possibility to disentangle compressional (EoS) 
and MDI effects on the basis of the present experimental results. It is worth to recall (see 
section I) that the theoretical predictions might also depend sensitively on other ingredi- 
ents in the model such as the in-medium nucleon-nucleon cross section (ann) which is not 
investigated in the present work. The possibility to constrain the in-medium cJnn has been 
discussed in previous works • We start with a very short description of the model and 

recall some of the typical features in the IQMD version. Then we discuss the influences of 
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the detector biases on the extraction of the in-plane anisotropy ratio. Finally, we confront 
the data to the model calculations for different parametrisations of the nuclear EoS. 



A. The IQMD model 



The Quantum Molecular Dynamics (QMD) model |T3| , p!7| , |i3| , |55| is a n-body theory which 
calculates the time evolution of intermediate energy heavy ion reactions on an event by event 
basis. Nucleons are represented as gaussian wave functions whose centroids are propagated 
according to the classical equations of motion and their width parameter L is kept fixed. 
These wave packets are submitted to two and three body local Skyrme forces plus Yukawa 
and long range Coulomb potentials. An additional term describing the momentum depen- 
dence of the nuclear interaction (MDI) is included. The latter, which is computed from 
the momentum dependence of the optical p-nucleus potential, produces an additional trans- 
verse deflection of nucleons early in the reaction. The parameters of the density dependent 
Skyrme and MDI potentials are adjusted to reproduce the ground state properties of the 
infinite nuclear matter and to fix the choice of the incompressibility K of the EoS. The 
numerical simulation of a complete event includes the initialisation procedure of projectile 
and target (see below), the propagation of nucleons through the previous potentials and the 
nucleon-nucleon interactions via a stochastic scattering term. For each binary collision, the 
phase space distributions in the final stage of the scattering partners are checked to obey 
the Pauli principle, otherwise the collision is blocked. At the final stage of the reaction 
{t = 200/m/c), composite particles are formed if the centroid distances are lower than 3/m. 

The IQMD (Isospin Quantum Molecular Dynamics) is a QMD version which takes 
into account isospin degrees of freedom for the cross sections and the Coulomb interaction. 
This version uses a gaussian width fixed to 4L = 8.66/m^ for Au nuclei. The initialisation 
procedure consists in distributing randomly into a spherical phase space the position r and 
the momentum p of the gaussian centroids with r < R and p < pp, where R = 1.12^4^/^ (A 
being the total nuclei mass number) and pp = 268MeV/c (Fermi momentum). 
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For the present work, the choice of this version, among other different reahsations of the 
QMD model, is justified by the fact that IQMD is particularly suited for the analysis of 
nuclear matter collective effects ||56|-pT|. Our study is however relevant in a more general 
level since the major deviations about nuclear matter flow between the different QMD ver- 
sions are now known to stem exclusively from the initialisation procedures and the gaussian 



widths [22 



B. Details about the calculations 



IQMD events were computed [^] for three different parametrisations of the nuclear EoS : 



i) a hard EoS {K = 380 A MeV) without MDI (H), ii) a hard EoS {K = 380 A MeV) with MDI 
(HM) and iii) a soft EoS (K = 200AMeV) with MDI (SM). The calculations were performed 
for the full range of impact parameters from to 14/m. The number of simulated events was 
about 300-1000 per fm interval. It is important to mention that the theoretical events were 
triggered using the same centrality selection as for the data. In addition they were filtered 
through the acceptance of the FOPI/Phase-I detector including geometrical cuts and energy 
thresholds. The calculations were restricted to charged particles with 1 < Z < 12. Here, 
the azimuthal distributions were determined with respect to the true reaction plane (known 
in the model). 

Before confronting the data to model calculations, let us make some comments on how 
the extraction of the anisotropy ratio R is influenced by the detector cuts. For this purpose, 
we used the IQMD (for both HM and SM parametrisations) events to compare the anisotropy 
ratios R before and after passing through the apparatus acceptance. The results are shown in 
Fig. ^ for the PM4 multiplicity bin, in terms of the ratio of unfiltered to filtered anisotropics, 
as a function of the cm. polar angle. As clearly evidenced by this figure, the distorsions 
caused by acceptance effects appear for very small 6c.m. and for 6c. m. > 65°. The values 
of R at small (large) 6c.m. angles are mainly affected by the Giab = 1-2° (©lab = 30°) 
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cut. Within ~ 4° and ^ 65°, the ratio R is only very httle affected (less than 5% as 
one can see from Fig. |^) by the detector apparatus. It is worth to note that the effect 
is very similar for both HM and SM calculations. Similar observations were also made 
for the H parametrisation. In the following, the discussion will be restricted to the region 
4° < 0c.m. < 65° where a realistic confrontation of data to model predictions can be done 
with a high level of confidence, independently on how well the experimental filter simulates 
the different apparatus biases. 

Using the same IQMD simulations (with a HM parametrisation), we have also investi- 



gated [0 the reliability of the method used in the previous section to correct the experi- 
mental anisotropy ratios for fluctuations of the reconstructed reaction plane. The validity 
of this method depends on how well the < cos(A$r) > and < cos^(A$r) > quantities 
are estimated 0. The latter have been calculated according to the method used for the 
data, i.e. by subdividing randomly each event into two and calculating the half difference 
in azimuth (A<I>r) between the reaction planes extracted from the two sub-events. They are 
compared in Tab. |T| to the values obtained using the difference in azimuth ($r — $true) 
between the reconstructed reaction plane and the true one (known in the model). As one can 
see the two procedures lead to very similar results, except for the PM2 and PM3 multiplicity 
bins where differences of about 5% are found. This is due to the relatively small number 
of detected particles per event in such rather peripheral collisions. We have examined the 
influence of this effect on the extraction of the anisotropy ratio. We found that even in 
the least favourable case PM2, the value of the anisotropy ratio {R =3.34) obtained within 
~ 4° and ~ 65° by applying the procedure used for the data to the IQMD simulated events, 
remains close to the one (/2=3.09) extracted from the azimuthal distributions relative to the 
true reaction plane. 
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C. Data versus IQMD predictions 



Tab. 1^ shows the mean multiphcities of different charged ejectiles in the PM4 bin for 
the data and for IQMD filtered calculations using a HM parametrisation. As one can see 
the model underpredicts by about a factor of 2.6 the experimental IMF multiplicities. The 
effect is even more important for Z = 2 particles (a factor of 4.3). It is found to be most 
pronounced for particles emerging from the participant region of the collision PB| , P5| , ^ . Sim- 
ilar deviations are observed for both SM and H parametrisations. This failure of the QMD 
model in reproducing the fragment multiplicities was also reported in other works P5|,p7|. 



It might be related p3 artificial equilibration of the excited nuclear matter which in 

reality breaks up into several fragments. It is also known that the clusterisation probability 
depends strongly on the gaussian wave packets width L used in the model Therefore, 



in order to perform a meaningful comparison between data and model predictions {i.e. with 
observables independent on the clusterisation process), one needs to reconstruct for both 
data and theoretical events coalescence invariant quantities. This is done, as described in 
section III, by including in the calculations of the anisotropy ratio all detected particles each 
being weighted by its charge. According to this procedure, the contribution of a given nu- 
cleon to the resulting observable is assumed to be independent of the fact that this nucleon 
is detected free or bound into a cluster. It has been shown in a recent work [^] that the 
coalescence process accounts for many of the observed features of the phase-space densities 
of light fragments emitted from the participant region, and in particular the increase in 
sideward flow with fragment mass. All the anisotropy ratios R reported in what follows are 
presented in terms of such coalescence invariant quantities. 

Figure ^ shows the ratio i? as a function of Oc.m. for semi-central events belonging to 
the PM4 multiplicity bin. Data are compared to the predictions of the IQMD model for the 
three parametrisations HM, SM and H. Data points are plotted with the overall systematic 
errors which have been estimated to 5%. As can be seen, the general trends observed in 



16 



the data are qualitatively well predicted by the model irrespective of the parametrisation 
used in the calculations. Indeed, theoretical R values reach a maximum around Qc.m. = 15°, 
mainly caused by the deflection of the projectile remnants, and tend toward 1 when Oc.m. 
approaches 90°. Now looking more quantitatively at the results, the following comments can 
be made. First, around Qc.m. = 15° the experimental anisotropy ratio is better reproduced 
by a hard EoS without MDI (H). Both SM and HM parametrisations, including MDI effects, 
overestimate the data by about 40% and 60%, respectively. This is disturbing as the MDI 
should offer a more realistic description of the reaction since the proton-nucleus scattering 
data call for a momentum dependent interaction. Thus, the previous observation would sug- 
gest that the MDI is not properly treated in the model. Nevertheless it is important to point 
out that other parameters, like the in-medium cTnn or the ones related to the initialisation 
procedure, could also influence the calculations [^. More detailled simulations performed 
with different sets of parameters should provide further information about this observation 
which remains difficult to be clearly interpreted at the present level of investigation. In 
particular, additionnal calculations made with a soft EoS without MDI would be helpfuU 
to draw conclusions about the effect of the MDI with respect to the stiffness of the EoS. It 
should be also noted that the use of coalescence invariant observables might be inappropriate 
at these small angles as the reaction mechanism behind the break-up and the evaporation 



of spectators is totally different from the freeze-out of a hot nucleon gas [p!7| , p8| , p5 |. 

On the other hand, the in-plane anisotropy carried by particles emitted from the hot and 
dense central region of the collision {i.e. at large Oc.m. angles) shows a clear sensitivity to 
the stiffness of the EoS and almost no dependence to the MDI. Data are better reproduced 
by both hard EoS calculations (H and HM). The SM parametrisation leads to anisotropy 
ratios lower by about 15% than those predicted by the H and HM versions. This observation 
agrees qualitatively with the results of previous analysis done with other kind of in-plane 
flow observables |g3|,||,||,|| | . 



Note that, by performing a similar analysis for the lower incident energies (100, 150, 250 
and 400^4 MeV), we have obtained comparable trends but lower sensitivities, as expected in 
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ref. |2],|7 . 

Figure |^ shows the ratio i? as a function of the impact parameter. Experimental events 
were sorted here into 9 multiplicity bins (each 5 multiplicity units wide) from PMUL = 
20 to 65. Theoretical R values were extracted for the four large multiplicity bins PM2 
to PM5 because of the low statistics. The impact parameter h for the experimental data 
was determined with the IQMD model using the HM parametrisation (see Tab. |ip. The 
comparison between data and model predictions is shown for two 6c.m. regions : 4° < 
Gem. < 25° (Fig. 0.a) and 45° < 0c.m. < 65° (Fig. |7|.b). The latter regions were taken 
from Fig. ^ which shows a sensitivity of R to the MDl and to the stiffness of the nuclear 
EoS for small and large 0c.m., respectively. In the low polar angular domain (Fig. 0.a), 
the H parametrisation underestimates the experimental anisotropics, in particular at impact 
parameters above 5/m; while both calculations including MDl (HM and SM) overpredict the 
data significantly. At large polar angles (Fig. 0.b), all three parametrisations underestimate 
the experimental results. This disagreement is, however, less pronounced with H and HM 
calculations. 

The comparisons in Fig. |^ and Fig. ^ show that the general trends observed in the data 
are qualitatively predicted by the model irrespective of the parametrisation used in the 
calculations, but, on a more quantitative basis none of these three model versions (H, HM 
and SM) can consistently reproduce all of the experimental results. To understand these 
discrepancies, one needs more detailed theoretical investigations including explorations of 
the influence of other physical ingredients such as the in-medium nucleon-nucleon cross 
section and the sensitivity of the theoretical predictions to the input parameters (such as 
the gaussian width) used for the initialisation. 

As illustrated in Fig. |^, the present data can potentially serve as a testing ground for 
the theory and should allow for disentangling compressional effects and those related to the 
momentum dependence of the nuclear force. Indeed in peripheral collisions, the anisotropy 
observable R exhibits an almost pure sensitivity to MDl at forward angles (Fig. |^.a) where 
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spectator matter is predominant; while at large 9cm angles {i.e. participant region) R is 
purely sensitive to the stiffness of the nuclear equation of state (Fig. 0.b). This can be under- 
stood as follows. For central collisions, the overlap between nuclei produces a highly dense 
nuclear matter volume where the escaping particles are strongly representative of the degree 
of compression reached during the reaction. In addition, the high number of binary collisions 
between nucleons lowers their initial relative momentum and consequently the repulsive ef- 
fect of the MDI. This is best seen by looking at the ratio R at the more central collisions for 
large cm. polar angles (Fig. 0.b), i.e. by avoiding the contribution of the projectile rem- 
nants which remain influenced by the MDI since they feel a relative small compression as 
compared to the participant nuclear matter. With increasing impact parameters the overlap 
between nuclei {i.e. compressional effects) and the number of nucleon-nucleon collisions de- 
crease. Hence, in contrast to central collisions, the ejectiles feel a higher dependence to the 
MDI and a smaller dependence to the stiffness of the EoS. In this regard, it is worth noticing 
that we observe the most pronounced sensitivity to the MDI by looking at the R ratio for 
small Qc.ra. angles in peripheral collisions (Fig. |^.a) where particles are almost exclusively 
deflected by the initial contact of the colliding nuclei since compressional effects are quite 
low. 

Although a definitive statement about the incompressibility of the nuclear EoS remains 
still elusive at the present stage, such comparisons have, nevertheless, the merit to illustrate 
the possibilities opened with the availability of the highly accurate data measured in the 
present experiment as stringent test for microscopic transport models, and in particular, the 
ability to unravel compressional and MDI effects. 

VI. SUMMARY AND CONCLUSION 

In this article, we reported the first detailled measurements of the in-plane azimuthal 
anisotropy ratio R of light particles and IMF emitted in Au+Au collisions at an incident 
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energy of 600 A MeV. Corrections for the finite resolution of the reaction plane reconstruction 
were applied. The study of the ratio R was carried-out by scanning the phase space in terms 
of cm. polar angle portions. This was found to be better adapted to the present detector 
acceptance as compared to rapidity selections. The ©c.m. dependence of the anisotropy ratio 
R exhibits a well defined maximum, located at forward polar angles (Oc.m. — 15°), which 
reflects mainly the bounce-off of the spectator remnants. We found that the larger the size of 
the detected particle, the stronger the magnitude of the azimuthal anisotropy and, therefore, 
the momentum space alignment with respect to the reaction plane. 

Data were compared to the predictions of the IQMD model using three different 
parametrisations of the nuclear EoS. The comparison was performed in terms of the ra- 
tio R determined as a coalescence invariant quantity, in order to avoid the problem of the 
lack of clusterisation in the IQMD model. It has been restricted to the c.m. polar angle do- 
main 4° < ©c.m. < 65° where detector acceptance effects on R were estimated to be less than 
5%. We found that the general trends observed in the data are quahtatively predicted by 
the model irrespective of the parametrisation used in the calculations; but, on a more quan- 
titative basis, none of these three model parametrisations (H, HM and SM) can consistently 
reproduce all of the experimental results. We show that, although a definitive statement in 
favour of one of the three parametrisations cannot be made at present, such comparisons 
have, nevertheless, the merit to illustrate the possibilities opened with the availability of the 
highly accurate data measured in the present experiment, over a broad range of impact pa- 
rameters, as providing a stringent test for microscopic transport models. In particular, the 
ability to unravel compressional and MDI effects was clearly demonstrated : in the spectator 
region, the anisotropy observable R exhibits an almost pure sensitivity to MDI, while in the 
participant region R is purely sensitive to the stiffness of the nuclear equation of state. 

Further detailed comparisons to microscopic transport model predictions, including ex- 
plorations of in-medium effects on the nucleon-nucleon cross section, should help to provide 
deeper insight into the properties of hot and dense nuclear matter and the underlying equa- 
tion of state. 
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TABLES 

TABLE I. Accuracy on the determination of the reaction plane for events measured under the 
PM2 to PM5 multiplicity conditions. (t(A<I>) is extracted from a gaussian fit to the A$ distributions 
(see text). The limits of the different event multiplicity classes are given for both data (PMUL) 
and IQMD (PMULjQ]y|-j-)). For each event class, the mean impact parameter (< b >) and the 
corresponding r.m.s. deviation extracted from the IQMD model are given. 



Multiplicity bin 


PM2 


PM3 


PM4 


PM5 


PMIT 


16-31 


31- If) 


16-G2 


> 62 


PMULiQMD 


19-38 


38-57 


57-76 


> 76 


<b> {fm) 


8.9 ±0.9 


6.7 ± 1.0 


4.0 ± 1.4 


3.0 ± 1.1 



(t(A$) {deg.) 34.5 23.6 24.7 30.0 



TABLE IL In-plane anisotropy ratio R before and after correction for fluctuations of the re- 
constructed reaction planes. The data are shown for the 40° < Gc.m. < 50° cut. Errors correspond 
to systematical uncertaincies. Statistical errors are lower than 1%. 



Multiplicity bin 


PM2 


PM3 


PM4 


PM5 


R not corrected 


2.32 ± 0.23 


3.36 ± 0.24 


3.02 ±0.15 


2.51 ±0.13 


R corrected 


2.86 ± 0.29 


3.94 ±0.28 


3.55 ±0.18 


2.96 ±0.15 
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TABLE III. Values of < cos($r - $true) >, < cosA$r >, < cos^($r - $true) > and 
< cos^ A$R > obtained from filtered IQMD calculations. See text for more details. 



Multiplicity bin 


PM2 


PM3 


PM4 


PM5 


< C0S($R - $true) > 


0.93 


0.97 


0.92 


0.88 


< COS A$R > 


0.87 


0.94 


0.92 


0.88 


< COs2($R - $true) > 


0.90 


0.94 


0.89 


0.85 



<cos2A$R> 0.85 0.90 0.87 0.81 



TABLE IV. Mean multiplicities of different types of ejectiles emitted in the forward cm. hemi- 
sphere for PM4 events class. Data are compared to different filtered IQMD model simulations. 
Errors represent statistical uncertaincies including the r.m.s deviations of the distributions. 





Z = \ 


Z = 2 


Z = 3 


3 < Z < 12 


DATA 


44.03 ± 0.05 


10.06 ± 0.02 


1.72 ± 0.01 


2.61 ± 0.01 


IQMD/HM 


61.89 ±0.27 


2.35 ± 0.02 


0.45 ±0.01 


0.99 ±0.01 
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FIGURES 

FIG. 1. Illustration of the FOPI/Phase-I acceptance in the {p[^\y^^'') plane for Au+Au colli- 
sions at 600 A MeV (see text for more details) . 

FIG. 2. Azimuthal distributions, with respect of the reaction plane, including the contribution 
of all detected charged particles, each being weighted by its charge. The data are shown for the 
semi-central PM4 event class, for different 0c.m. bins reported on the plots. Data points are larger 
than the corresponding statistical uncertainties. The curves are the results of the fit described in 
the text. 

FIG. 3. In-plane anisotropy ratio R as a function of the polar angle Qc.m., for different types 
of ejectiles. The data are shown for the PM4 event class. The error bars are of statistical origin. 
The dotted curves represents a gaussian fit to the correlation between R and logio{Qc.m )- 

FIG. 4. Maximum in-plane anisotropy i?™*^ (top-panel) and the corresponding polar angle 
®™n? (mid-panel) as a function of the charge of the detected fragment. The lower panel shows the 
Z dependence of the mean multiplicity in the forward cm. hemisphere. The data are shown for 
the PM4 event class. Error bars include statistical errors only. The curves are just to guide the 
eye. In the lower panel, data points are larger than the corresponding statistical uncertainties. 

FIG. 5. In-plane anisotropy R obtained from unfiltered IQMD simulations divided by the one 
corresponding to filtered calculations, as a function of the cm. polar angle Gem.- The calculations 
are shown for Au(600AMeV)-|-Au collisions preselected with the PM4 multiplicity cut. The full 
(dotted) line corresponds to a HM (SM) parametrisation of the nuclear EoS. The horizontal dashed 
lines indicate i? = 1 ± 0.05 values. Statistical uncertainties in the IQMD values are lower than 2%. 
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FIG. 6. In-plane anisotropy ratio R as a function of the cm. polar angle for semi-central 
(PM4) collisions. Data (circles) are compared to the predictions of the IQMD model with the 
paramctrisations HM (solid curves) , SM (dotted curves) and H (dashed curves) . Theoretical values 
are calculated with respect to the orientation of the true reaction plane after application of the 
FOPI detector filter. Experimental values are corrected for the unaccuracy on the determination 
of the reaction plane (see text). The shaded zone corresponds to the polar cm. angular domain 
where the data suffer significant distorsions due to apparatus effects. Error bars on data points 
correspond to systematical uncertainties estimated to 5%. Statistical errors on theoretical points 
are lower than 2%. 

FIG. 7. In-plane anisotropy ratio R extracted in the cm polar angular domains 
4° < Gc.m. < 25° (a) and 45° < 6c.m. < 65° (b) as a function of the impact parameter (6) . Data 
are compared to the predictions of filtered IQMD model calculations. Error bars on data points 
correspond to systematical uncertainties estimated between 5% and 10% depending on the collision 
centrality. The impact parameter is deduced from IQMD calculations using a HM parametrisation. 
Statistical errors on theoretical points are lower than 2%. 
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